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SECTION 1.0 
INTRODUCTION 


The GEOS-3 spacecraft, launched April 10, 1975, was the first orbiting 
vehicle to carry a coherent C-Band transponder. The previous spacecraft in the 
GEOS series carried non-coherent C-Band transponders, which were not suitable 
for two-way Doppler range rate measurement. 

Approximately one-third of the C-Band radars available to support the 
GEOS-3 mission have coherent capability, either the Coherent Signal Processor 
(CSP) or Velocity Extraction Subsystem (VESS). These systems produce a digital 
measurement of target range-rate along the slant range vector. The measure- 
ment system is implemented as an additional Type II servo, modeling R and R and 
driven by error signals extracted from return signal frequency shift compared to 
the local oscillator, with suitable provisions for resolving harmonic ambiguity. 

Prior to the launch of GEOS-3, these systems were able to produce range 
rate data on orbiting vehicles only when the radar was operated in the skin track 
mode. This data was characterized by severe noise and signal strength problems. 
The advent of a spacecraft-borne coherent transponder allowed the examination 
of precise C-Band range rate data in the context of a trajectory highly 
constrained by orbit dynamics. The general quality of this data is illustrated by 
Figure I- 1, showing range rate residuals for a single GEOS-C pass taken by the 
AN/FPQ-6 radar at Wallops Island, Virginia. The residuals are about an orbit fit 
to the range and angle data taken by the radar. This data shows low noise (0.65 
cm/sec RMS) and little systematic error after correction in the data reduction 
process for known timing and measurement bias errors. The identification and 
elimination or reduction of these errors is the subject of this report* 






The availability of two independent data sets simultaneously from the 
same instrument suggested that calibration efforts might be enhanced by the 
ability to compare them in the same measurement domain. The idea of 
integrating range rate data into pseudo range measurements is not new; a similar 
study was undertaken in support of GEOS-2 by Wells and Guard (I). Their results 
were severely limited by the fact that only range rate data from the skin 
tracking mode was available. Using this work as a foundation, we have extended 
the technique using the more precise range rate data available from the GEOS-3 
coherent transponder. 

Results from the integrated range evaluation led to the development of 
additional corrections to the conventional C-Band ranges, leading to a reduction 
both in noise and systematic error. These techniques include: 

• Filtering of data to reduce relatively high radar recording rates (10 
or 20 pps) to convenient rates for orbital data reduction (usually I 
point per 6 seconds). 

• A radar AGC based correction for beacon delay variations with 
signal strength. 

The former reduced typical RMS range noise figures for the Wallops Island FPQ-6 
from one meter or more to the 25-50 cm range, while the latter corrected for a 
systematic error whose dynamic range through a typical pass can approach 2 
meters. 


SECTION 2.0 

PROCESSING TECHNIQUES 


/ This section details the processing techniques used to produce and 
analyze the integrated Doppler range data, including the preprocessing correc- 
tions necessary to the conventional* ranges produced by the radar, which define 
the reference orbit for evaluation of integrated range results. The overall data 
processing organization is shown in Figure 2-1. 

All of the data presented in this report were obtained by the AN/FPQ-6 
radar at Wallops Island, Virginia. This radar is equipped with the Coherent Signal 
Processor pulse Doppler modification and a DIRAM range machine. Tracking data 
is recorded digitally in real-time under control of the radar’s RCA 4101 
computer. The remainder of the data processing activities were carried out at 
the Wallops Flight Center Information Processing Laboratory, on either the 
Honeywell 625 or Honeywell 635 medium scale computer systems. 

The initial step in the processing chain was the reformatting of the 
binary site tapes into a form compatible with the 625/635 computer system. This 
step is undertaken by the PASS I program* The output of this program is then 
input into the PASS2 C-Band preprocessor. It is within this program that the 
filtering, integration, and other calibration techniques discussed in this report 
were implemented. The output of the PASS2 program, consisting of the final 
reduced and corrected data, was then used in the orbit determination program 
GEQDYN. Resulting residual information was passed to a residual analysis 
program GEORGE, which was modified to incorporate appropriate integrated 
range error models. In parallel ta this sequence, the ORAN error analysis 
program was used to predict the effects of the various modeled and unmodeled 
errors on the batch least squares orbital solution. 


In the case of the Wallops Island AN/FPQ-6 these are often referred to as 
DIRAM ranges, after the Digital RAnge Machine, which is the range tracker on 
that radar. 


R, R, 
TIME 


PASS 1 
REFORMAT- 
TER 


PASS 2 

PREPROCESSOR 


INTEGRATE RATES, 
REFORMAT, TRANSIT 
TIME 




GEODYN 

ORBITAL 

SOLUTION 



1 GEORGE 
RESIDUAL 
ANALYSIS 



DATA 

REDUCTION 

RESULTS 




FIT ORBIT TO RANGE 
AND ANGLE DATA; 
GENERATE CSP 
RANGE RESIDUALS 


PLOT RESIDUALS, 
FIT ERROR MODELS 


ORAN 

ERROR 

ANALYSIS 


ERROR 

ANALYSIS 

RESULTS 


SIMULATE ORBITAL 
SOLUTION WITH 
MODELED AND UN- 
MODELED ERROR 
PROPAGATION 


FIGURE 2-1. INTEGRATED RANGE PROCESSING AND ANALYSIS 








The relevant portions of each of tht«e processing steps are discussed 
Individually in the following sections. 


2. 1 INTEGRATED DOPPLER RANGES 

The range rate measurements at the radar sampling rate of 10 or 20 pps 
are numerically integrated with a simple trapezoid rule integration scheme 
formulated as, 


R i + P 'M 

R - R + -i — 4t 

I I" I ^ 


This scheme is graphically illustrated in Figure 2-2. Despite Its 
simplicity, this algorithm has several advantages. In particular, it 

• does not require equally spaced samples in time, 

• is fast and efficient computationally, and 

• has error dependent on the second derivative af the integrand (here 
IT), i.e.j the non-linearity of the function to be integrated between 
consecutive data points. 

In general, this scheme has been found adequate for the study at hand. 
The function (range rate) is sampled frequently enough (10 or 20 pps) in relation 
to the magnitude of the dominant error term (R) to limit integrator error to a 
few cm in magnitude. A more complete analysis of the integration technique is 
included in this report as Section 2.4.2. 
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The determination of the constant of integration, which in this case 
amounts to a "zero-set" or bias correction for the integrated ranges, is initially 
accomplished by adjusting them in a batch sense to the ensemble of DIRAM 
ranges. The integration process is begun using the DIRAM range associated with 
the first valid range rate measurement. As the process progresses, the sum of 
the differences between the integrated Doppler range and DIRAM range at each 
point is kept. At the end of the pass the ensemble of integrated ranges is 
corrected by the mean difference, i.e., 



where 

R = ith C-Band DIRAM Range 

i 

R^ = ith integrated Doppler Range 

i 

In effect this introduces into the integrated range data set a bias identical to 
that in the DIRAM range machine. 

The integrated ranges are sampled at every 60th point for an output data 
rate of once every 6 seconds. This is a convenient rate for data in orbital 
solutions and corresponds to that chosen for reduction of normal range, range 
rate, and angle data. Because of the low noise on the data, there is minimum 
benefit to be gained by filtering the integrated ranges. 

It should be noted that prior to the integration process, range rate data is 
corrected for bias using an estimate determined by a pre- and post-mission "one 
line up - one line down" calibration procedure, which is further discussed in 
Section 3.0. Additionally, measurement time tags are corrected for propagation 
delay and the data is reformatted suitably for input to an orbit determination 
program. 
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2.2 DIRAM RANGE CORRECTIONS 

The principal preprocessing corrections now made to the DIRAM ranges, 
in the order of their development and implementation, are 

• Bias adjustment based on pre- and post-mission ranging to a 
surveyed target. 

• 60 to I reduction of sampling rate through application < r a 
polynomial midpoint filter. 

• Correction for within pass beacon delay variation based on radar 
receiver AGC. 

The first of these, bias correction based on pre- and post-mission ranging 
to a surveyed target, is a conventional technique long used in a great many 
ranging systems. The Wallops Island FPQ-6, with suitable pre-mission warmup 
time, does not drift significantly in range bias during the mission interval. The 
remaining corrections were first implemented during the GEOS-3 study, and are 
discussed individually. 


Filterinc 


C-Band radars typically record data at either 10 or 20 pps, a considerably 
greater rate than is necessary or useful for orbit determination and geodetic 
studies. Historically this data has been reduced by sampling; i.e., the often used 
once-per -six-second data rate was achieved by selecting every 60th or 120th 
point and ignoring the rest. It was felt that a relatively simple filtering 
technique might considerably reduce the noise of the resultant data. The intent 
was not to derive an optimal or near optimal filter or smoother for the range 
data, but rather to find a simple and computationally efficient method to stop 
"wasting" most of the information recorded. 
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The solution implemented was to apply a quadratic polynomial midpoint 
filter to each set of 60 or 120 points. Serial correlation effects are minimal 
since the filter is not applied recursively but rather to each successive six second 
interval of data. In effect, it is simply a one higher order technique than 
averaging each consecutive set of 60 points. Moreover, the filter is implemented 
as a set of 60 weighting coefficients, which does not add significantly to the 
computational burden of preprocessing the large amounts of C-Band data 
processed at Wallops Flight Center, Figures 2-3 and 2*4 show the results of 
filtering on a typical range data pass. 


2.2.2 Beacon Delay Variation 

Both the coherent and the non-coherent C-Band transponders on the 
GEOS-3 spacecraft were the subject of detailed prelaunch calibration and test 
procedures as reported by SGLSER (2). Figure 2-5, excerpted from his report 
shows the delay variation of the coherent transponder with received signal 
strength. In the case of the AN/FPQ-6 at Wallops Island, with the normal GEOS- 
3 setup, received signal strength at the spacecraft transponder can vary over the 
range 4Q to -30 dBM within a single pass as the range varies from almost 3,000 
Km at 5 elevation to less than 1,000 Km at high elevation angles, As the figure 
illustrates, this causes a variation of approximately 6ns in two-way propagation 
delay or a change of over a meter in the radar measured range. This significant 
systematic effect was uncorrected in the past due largely to the difficulty of 
extracting transponder received signal strength information from the spacecraft 
TM stream and suitably correlating it with the C-Band data processing effort. 
Nevertheless, proper evaluation of the extremely low noise integrated Doppler 
ranges could not be completed with a systematic error of this magnitude in the 
ranges determining the reference orbit. Accordingly, a scheme was developed 
whereby the digitized AGC level at the tracking radar could, with suitable 
calibrations, be converted to radar received signal to noise ratio, and ultimately, 
through the radar equation and knowledge of the radar's and spacecraft's 
transmitted power, to spacecraft transponder received signal strength and 
finally, beacon delay variation. This process is illustrated in Figure 2-6. 
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FIGURE 2-6. RADAR AGC BEACON DELAY VARIATION CORRECTION PROCESS 





Initially the radar performs an AGC-stepcal procedure which functionally 
relates the radar recorded AGC level to the radar received signal-to-noise ratio. 
The radar AGC from each point during tracking can then be converted, through a 
table-lookup linear interpolation process, to radar received SNR. This in turn 
can be converted to transponder received signal strength through a K-factor 
computation (Figure 2-7) making use of the radar equation and knowledge of 
radar and transponder operating parameters. 


Finally, beacon received signal strength is converted to relative delay by 
another table-lookup linear interpolation based upon the pre-mission test results 
reported by SELSER. The scheme is implemented with a relative beacon delay of 
0 (i.e., nominal delay of 379.48 m) at the -30 dB beacon received signal strength 
level. 

2.3 ORBIT DETERMINATION 

Reference orbits for all of the analyses reported here are single station, 
single pass trajectories fit to the DIRAM range and angle data from each radar 
pass. The orbit determination program used was the GEODYN program (3). 

Single station single pass orbital solutions almost always fit very well, 
the orbit largely accomodating most systematic errors, particularly station 
position errors and absolute timing errors which are near unobservable in data 
from a single station. This method for determining the reference trajectory for 
integrated range analyses was chosen specifically for these reasons. The low 
noise integrated Doppler ranges would be very sensitive to almost all such 
unmodeled errors. By using a reference trajectory determined by an identically 
colocated station using a common timing standard and position, timing and 
gravity model errors are essentially eliminated. Furthermore, tropospheric 
refraction correction error, identical for both sets of ranges, is largely 
accomodated by the orbit. The remaining systematic errors observable in the 
integrated range residuals about such a reference orbit are listed in Table 2-1. 



FIGURE 2-7. K-FACTOR COMPUTATION 
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TABLE 2-1. 

INTEGRATED RANGE ERROR SOURCES 
ERRORS VISIBLE IN INTEGRATED RANGE RESIDUALS 

I. SYSTEMATIC ERRORS IN RANGE AND ANGLE DATA USED TO DETERMINE 
THE ORBIT 

• RANGE BIAS 

• RANGE TRACKER DYNAMIC LAG 

• IONOSPHERIC REFRACTION 

• TRANSPONDER DELAY VARIATIONS WITH SIGNAL STRENGTH 

II. SYSTEMATIC ERRORS IN THE COHERENT SIGNAL PROCESSOR 

• MEASUREMENT BIAS 

• SERVO LAG 

• IONOSPHERIC REFRACTION 

• RELATIVE TIMING ERROR 

COMMON ERRORS LARGELY ACCOMMODATED BY THE ORBIT 

• TROPOSPHERIC REFRACTION CORRECTION ERROR 

• ABSOLUTE TIMING ERROR 

• STATION POSITION ERROR 




Category I are orbit errors caused by systematic errors in the orbit- 
determining ranges but not present in the integrated range data. 

' ) ’ ; %/' ■■ V 

Category li are data errors in the CSP integrated ranges not correspond- 
ing to errors in the DIRAM ranges. Ionospheric refraction appears in both 
categories since it has equal but opposite effect on the group-measurement 
DIRAM ranges and the phase-measurement integrated Doppler ranges. 

2.4 ERROR ANALYSES 

CSP integrated range error analyses were undertaken for three distinct 

sources: 

• radar systematic errors, 

• computational errors: integrator error, 

• uncorrected propagation effects: ionospheric refraction, tropo- 

spheric refraction correction error. 

These analyses were conducted both deductively through model fitting to 
actual residuals and inductively using the ORAN program (4) to simulate the 
effects of unmodelled errors in a Bayesian Least Squares orbital solution. 


Integrated Range Error Model 


After extensive investigation, the most reasonable form for modelable 
error in the residuals of the integrated range data about the orbit determined by 
the conventional data was found to be: 



4R i = a o + a l R i + a 2 (t r’o ) + a 3 R i 
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where 


4Rj a ith integrated range residual 

a Q a offset or bias; error in the constant of integration 

a, s relative timing error between the CSP and DIRAM systems 



Each of the four error terms represent systematic errors known to exist and 
quantitatively modelable. This model was fit in a batch least squares sense 
(Appendix B) to integrated range residuals generated in the orbit determination 
program GEODYN, 



2.4.2 Integrator Error 


The integrator for CSP ranges is a simple trapezoid rule formulation 
shown in Figure 2-2. It takes its name from the fact that the area of the shaded 

. *, t. 

rectangle is the area of the trapezoid (t a R,_ j, Rj, t.). Aside from its ease in 
implementation, this integrator has the following characteristics: 
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• It is precisely accurate if the integrand Is linear between adjacent 
sample points; thus its first order error is a function of the second 
derivative of the integrand. Since the second derivative of range 

•H ^ 

rate (R or "jerk") typically does not exceed .30 m/sec J and the 
j integrator is exercised at the 10 pps radar sample rate (step size h 

i - ' . • 

j ~ . I sec) the local error written (5) as 

t ' 

j . . ... ■■ • . ■ • 

! - -L h 3 f" (£) 


i is not expected to exceed .0025 cm. and the global error 

I t. - t ? 

- J — - fr f" (e) 

12 

j. . 

should not exceed 15 cm., assuming that the maximum R value 
prevails throughout the pass. 

!■• ■ • ■ ■ ' . . ' • • ... ■ . . , 
■■■■■■ • •• • ■ • ■ 

• It does not require equally spaced data points; it can integrate over 
editted points and, moreover, can be used on the series of unequally 

j 

spaced (transit-time corrected) measurements generated by the 
radar. 

i • • ■■•••; •...■• • • ■ . 

i /■;. 

! .. • 

It is essential to the integration process that the sample times be 

corrected for propagation delay (transit time) effects before the intersample *it 
is computed; this intersample period is not constant as the propagation delay 
changes with range. Accordingly an integrator able to function with unequally 
spaced data is a necessity, unless a measurement interpolation scheme is 
implemented. The following clarifies the necessity of this correction. 
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The numerical integrator functions upon an ordered set of (t., R.) values 
as sampled by the radar. The measurement times are at even .1 sec intervals of 
radar (ground received) time. If we assume that these measurements are 
"instantaneous" then we must alter the time associated with each measurement 
to its true effective time at the satellite, 


where 


Tj = true effective time of ith measurement 
t. = ground time of ith measurement 

Rj = range at ith measurement 


Whether the times associated with the measurements are altered before 
the integration process (while associated with a rate measurement) or after the 
integration process (while associated with a range measurement) is irrelevant in 
principle since the integral of ground-received rate data is ground-received range 
data, and conversely the integral of satellite-timed rate data is satellite-timed 
range data. This theoretical indifference to transit time considerations rests 
upon the notion of perfect theoretical integration. In fact the numerical 
integrator we use makes finite approximations of area at each step and to do so 
accurately must have a bin width (4T) equally as accurate as the bin height (R). 
Even so, it too is indifferent to a time bias, i.e., it calculates 
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and if all times are in error by a bias 


*i - T i ♦ « 


we still have 


= f i ' Vi 


= (T. + €) - (T._, + e) 


dT 


The fundamental point here is that transit time is not a bias and does not 
fit this model. Even for two points .1 sec apart, the true ’'width*' of the "bins" 
differs significantly. Mathematically, 


R. 


T i = *i - 


R. 


V. = Vi 


i-l 
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T i - T i-I 


= (t. - t. ,) 

i i- 1 


R i- R i-I 
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which simply says that our bin widths differ from .1 sec by an amount 
proportional to the change in range over .1 sec. This is clearly a range rate 
dependent error; indeed, we can write 


AT 


AT 


vw - 


(R; - R;_ 


<*i ->!-!> 


I 

I 

l 



4T - (♦rVtX'-J* 

as (t. - tj j) — » 0 


The magnitude of this error is clearly illustrated in Figure 2-8, a typical 
GEOS-3 range rate plot with two additional scales, transit time rate, and actual 
intersample AT at the 10 pps ground received sampling rate. 


Returning to the integrator formulation, 
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R. +R 


i-l 
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we find that an error in 4T of produces a local error in the integrated 
range increment of 



R^ dt 


6 « 

c 

Rj 2 dt 

Figure 2-7 shows the local error — p- - for a 55° pass. The expected 
global error for the pass is the sum integral) of these local errors. Figure 2- 
10 shows the actual residual plot of Rev 4547 integrated with a constant .05 sec 
(20 pps) step size, not corrected for transit time. This error is independent of 
stepsize, since doubling the stepsize results in half as many steps, each with 
twice as much error. 

2.4.3 Propagation Effects 

2*4.3. 1 Tropospheric Refraction Error 

As previously stated in Section 2.3, tropospheric refraction effects 
should not produce significant errors in the integrated range residuals as 
processed for this study. Ranges are corrected for tropospheric effects with a 
model (6) incorporating local meteorological data and good to approximately 2- 
3%. The residual refraction error is a smooth function, primarily of elevation, 
which is almost entirely accommodated by the reference orbit fit to the 
conventional range data at elevation angles above 15°. Since the residual error 
in the integrated ranges is identical, it is not observable in residuals about the 
reference orbit. 
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2.A.3.2 Ionospheric Refraction Error 

Accurate quantitative modeling of ionospheric refraction error effects is 
extremely difficult and can require reasonably good estimates of solar flux 

■ . . . i _ 

values for the time of the data in question. In general, C-Band tracking data has 
not been routinely corrected for ionospheric effects by any investigators known 
to us. However, since the effect of the ionosphere will be equal but opposite on 
the two ranging data set types discussed here, and because of the very low noise 
of the integrated range residuals where the effects were expected to be visible, 
some quantitative modeling was attempted using a modified version of the BENT 
ionospheric model in the ORAN program. The results are discussed in Section 
3.0. 
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2.4.4 Range Bias Effects 

Analysis of the pass presented in detail in the following section, Rev 
9799 as tracked by the Wallops AN/FPQ-6, benefited from a reliable estimate of 

the total range bias in the DIRAM ranges used to compute the reference orbit. 

! - ’ 

This bias was recovered from a two-consecutive - revolution solution making use 
of the tracking data available from the same radar for Rev 9800. During the 
course of the study, however, many passes were examined which either had no 
such estimate or which were initially preprocessed incorrectly in such a manner 
as to introduce a large bias. In each of these cases a distinct "M" or "W” shape 
was observed in the integrated range residuals and where the bias was extreme 
(such as neglect of transponder delay, « 380m) in the ordinary DIRAM range 
residuals. 


An ORAN analysis determined that this shape is characteristic of the 
effect of a range bias on a single-station single pass orbital solution. The sense 
of the pattern (M or W) is determined by the sign of the bias. It is best explained 
as jthe geometrical resolution of a constant bias in slant range into varying 
proportions of aiong-track and radial error, as the station-satellite geometry 
changes during the pass. The pattern is produced by the nature of a least-squares 
orbital solution and its attempt to find the one trajectory that best accomodates 
this changing error. Figure 2-1 1 shows the effect of a 10m bias as predicted by 
the ORAN program for both a low and high elevation pass. As might be 
expected, the effect is most pronounced in high elevation passes. 

In fact this phenomenon had been noted before in cases where an 
extremely large bias had been introduced into conventional range data through 
error or oversight. Its significance was minimal because of the magnitude of the 
effect: as the figure shows, a 10m bias in a high elevation pass produces only a 
Im effect; in prior analysis of C-Band ranging data with an RMS noise level of 
lm or more such an effect was largely unresolvable. However, the availability of 
extremely low noise range measurement from a C-Band radar may significantly 





increase the importance of this phenomenon. Scaling from the figure, a Im bias 
produces an "M" or "W" shape on the order of IQ cm at high elevations; this is 
we|l below the potential resolution of the integrated Dbppler range data. With 
the availability' of accurate radar error parameter values and good quantitative 
estimates of ionospheric propagation effects, the possibility of recovering range 
bias estimates to less than a meter from data at a single station during a single 
pass could significantly enhance the value of the C-Band network in geodetic 
studies. 
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SECTION 3.0 
RESULTS 

i . ■„ ' . 
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I The results of the rather complex chain of processing, correction, error 

modeling, and analysis of integrated CSP ranges and the equally important 
reference orbit determining conventional DIRAM ranges is described sequentially 
| for a typical data set in Section 3.1. This somewhat iterative solution is also a 

; semi-chronological recreation of the development of the various techniques over 

j the course of the study. Representative results from other passes are presented 

!■:. in Section 3.2. 

I 

j Prior to any of the analyses presented here, rate data and the integrated 

j ■ ' ; 

ranges produced from It are corrected for a constant timing error relative to the 

conventional range data. This error was discovered and analyzed prior to the 
j work described in this report and has been described in detail by Krabill and 

Dempsey (7). It is related to the update rate of the local oscillator in the CSP 

system. The error was determined to be correctable as a timing bias equal to 

one-half the update interval, usually the radar's pulse repitition interval (PRI), 
and has become known as the error. It is not discussed further in this 
report,. 


3.1 REV 9799 



Figure 3-1 shows the results of applying the basic integrated Doppler 
range technique to Rev 9799 as tracked by the Wallops Island AN/FPQ-6 on 
March 2, 1977. No special corrections were made to the range rate data before 
or during the integration process other than the timing error previously 
discussed. The simultaneous DIRAM range data set is shown in Figure 3-2. It has 
been corrected for bias on the basis of pre- and post-mission ranging to a 
surveyed range target. Both data sets have been corrected for transit time and 
tropospheric refraction effects. The DIRAM ranges were reduced to the 10 per 
minute rate by the polynomial midpoint filter technique detailed in Section 2.2.1. 
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The principle trend in the uncorrected integrated range data is the 

I j ■ : I" i • - - 

monotonic decrease consistent with a constant bias (~ 2 cm/sec in this case) in 

the range rate measurement. It illustrates the sensitivity of the integrated 
ranges to systematic error and their ability to magnify such effects. The bias is 
a known error stemming from a discriminator curve offset. A hardware 
calibration technique was developed to estimate this bias, based on driving the 
CSP system up one spectral line and down one spectral line during pre- and post- 
mission calibration periods. This data is recorded on tape and can be reduced; to 
yield the true CSP zero, which is essentially the mean of the two e ->ual in 
magnitude but oppositely signed rate measurements. Figure 3-3 shows the data 
after this calibration correction to the range rate measurements is taken into 
account during the integration process. 

The residual effects were thought to be some combination of; 

• Servo dynamic lag 

• Ionospheric refraction effects 

• DIRAM range bias effect on the reference orbit 

• Residual timing bias in the CSP system (after PRI/2 correction). 

However, attempts to recover (Appendix B) consistent estimates of error 
model parameters (lag coefficient and timihg bias) from this and other data sets 
proved futile. The various error sources appeared to produce highly correlated 
effects that were difficult to separate in a least squares recovery, i.e., the 
solution was weak and the estimates highly correlated. Most disturbing of all, 
the characteristic residual patterns seemed not to show the predicted effect of 
servo lag in the conventional sense. In fact, the residual patterns and the lag 
coefficient estimates recovered! from them tended to be the reverse of what was 
expected, i.e., the servo appeared to "lead." This observation initiated a search 
for a systematic error source whose signature might be nearly opposite to that 
expected from CSP servo lag effects. This error source was ultimately identified 
in the DIRAM ranges defining the reference orbit. 
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3. 1 .2 Beacon Delay Variation 


As described in Section 2.2.2, pre-flight test data indicated that the 
GEOS- 3 coherent beacon delay might vary significantly with signal strength in 
nominal operating conditions. Analysis indicated that this error effect in the 
DIRAM ranges, if accomodated by the reference orbit, could produce an "error" 
pattern in the integrated range residuals quite similar to that expected from CSP 
servo dynamic lag effects, but of opposite sign. Accordingly, the technique of 
correcting range data for beacon delay variation was developed and applied. 

Figure 3-4 shows the range correction this technique produced for each 
of the DIRAM range measurements during the pass in question. This figure shows 
that, as calculated by the K-f actor method previously described, the beacon 
delay was approximately nominal 379.48 m at PCA for this pass. To either side 
of PCA, the beacon delay was longer due to decreased signal strength at longer 
ranges; the proper compensation being a negative range increment to "shorten” 
the recorded ranges - in this case by as much as two meters. The step-like 
effect prior to PCA represents toggling of the relatively coarse radar AGC least 
significant bit. The large excursion after PCA is probably representative of an 
antenna null or lobing effect. The overall bell shape could easily be accomodated 
by the reference orbit fit to this single pass from a single station. The excursion, 
however, should have been evident in the DIRAM range residuals. Figure 3-5 
shows the uncorrected DIRAM residuals and the residuals after correcting the 
data for beacon delay variation and fitting a new reference trajectory. The 
reduction of the group of "long" ranges in the area of the antenna null event is 
apparent. 

Figure 3-6 shows the same integrated range data as Figure 3-3, displayed 
however as residuals about the reference orbit determined by the beacon delay 
variation corrected DIRAM ranges. A "bucket" shaped pattern more consistent 
with the expected effects of CSP servo dynamic lag is visible. Least-squares 






METERS 


FIGURE 3-5. 

REV 9799 RANGE RESIDUALS 
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estimation of error model parameters yields the results shown in Figure 3-7. 
These results show considerable improvement over those that had been obtained 
earlier - essentially negligible residual timing and range rate bias errors, and a 
lag coefficient estimate K Q = 165, compared with the nominal K q of 225 for the 
servo in this bandwidth. The residua! data illustrates the extremely low noise (3 
cm RMS in range ) of this data. However, the data cannot be considered fully 
correct as a significant ionospheric refraction effect should still be present. We 
conclude that the multi-parameter error model can accomodate almost any 
smooth function through aliasing in one or another of its parameters. Although 
useful in the exploratory stages of this technique, it was concluded that 
consistent a priori estimates of system error parameters were necessary. 

i • i : ’ v : ■ . ■ ■: ■ . 

The nominal estimate of system relative timing error after the PRI/2 
correction is essentially zero (7)* The estimation of range rate bias with the one 
line up/one line down technique is thought to be quite good; however a 10% error 
in the estimation of a typical bias (1.5-2 cm) would significantly affect 
intjegrated range results. Similarly, the nominal servo lag coefficients K q , are 
expected to be accurate + 15%, but this also is too great a range for consistant 
data reduction of the integrated ranges. For example, Figure 3-8 shows the 
results for the subject pass when corrected, a priori, for 

• Residual timing bias: 0. 

• Range rate bias: one line up - one line down estimate of -2.04 

cm/sec. 

• Servo lag: K = 225. 

a 

This can be compared to Figure 3-9, the predicted effects of ionospheric 
refraction as observed in the corrected integrated range residuals about the 
reference orbit. The disagreement between these two figures could be explained 
almost entirely on the basis of uncertainty in the error model coefficients and 
difficulty in obtaining accurate quantitative estimates of ionospheric propagation 
effects. 
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3.2 OTHER CASES 


Several other typical results are shown in Figures 3-10 through 3-12. 


Figure 3-10 shows the integrated range residuals for Rev 202 after the 
least squares solution for the radar error parameters. Again, the values 
recovered are in reasonably good agreement with the nominal estimates of K = 

I • 

27, | R bias after calibration = 0., relative timing error after PRI/2 correction = 0. 
However, the effects of ionospheric refraction have been largely absorbed by 
these estimates. In particular, the relative timing error estimate of almost half 
a millisecond is much too large and probably due in large part to this effect. The 
remaining residual pattern shows a "W" shape consistent with a small bias in the 
DIRAM ranges used to compute the reference trajectory. The small tails at the 
extremes of the "W" are thought to be the effects of tropospheric refraction 
model errors which are not totally accommodated by the reference orbit at very 
low elevation angles. The effect here is limited to data between 10° and 15° 
elevation. 


Figure 3- 1 1 shows the integrated range residuals for the next consecutive 
revolution, Rev 203, corrected a priori for the nominal error parameters as 
indicated in the figure. The remaining pattern is characteristic of predicted 
effects of the ionosphere, which are not always symmetric due to spacecraft 
motion toward or away from the sun. Deviation of the pattern from the actual 
effects of the ionosphere is attributed to error in the a priori estimates of radar 
error parameters used to correct the data. 


Figure 3-12 shows the integrated range residuals from Rev 1 1704, once 
again after correction by least-squares recovery of error parameters. In this 
case the effect of the ionosphere apparently served to largely cancel servo lag 
effects; the recovered K q was essentially infinite indicating no lag, while the 
range rate timing and measurement bias estimates were more reasonable. 
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FIGURE 3-10. REV 202 INTEGRATED RANGE RESIDUALS 


Ji SOLVED PARAMETERS K_ = 32 

3 
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FIGURE 3-12. REV 11704 INTEGRATED RANGE RESIDUALS 
SOLVED PARAMETERS K a = « (no lag) 

R BIAS = .7 mm/sec 
At = —143 /me 



Other passes studied tended to yield similar results; the integrated 
ranges could always be reduced to a very low noise (1-3 cm RMS) data set by 
fitting an error model whose terms accommodated the effects of the ionosphere. 
Conversely, correction on the basis of a priori error term estimates yielded 
systematic error patterns grossly similar to those predicted by cursory 
ionospheric modeling. Uncertainty levels in both the ionosphere and the radar 
models indicate that the results are consistent. 



SECTION 4.0 
CONCLUSIONS 

The primary objective of the GEOS-3 spacecraft C-Band Experiment was 
defined in the GEOS-C Mission Plan as, 

"To better determine the absolute accuracy of instrumentation 
radar systems, develop refined methods of calibrating these systems, and 
improve the techniques employed in processing the associated data.” 

This study has directly addressed that objective and made significant contribu- 
tions to its achievement. 

The integrated Doppler range technique has proven itself a sensitive and 
useful tool for resolving systematic errors and aiding in the calibration of 
coherent C-Band radars. In addition to producing extremely low noise range 
data, the potential for resolving the systematic error pattern present in a single 
station single pass orbital solution using biased range data may lead to an 
entirely new technique for calibration of C-Band radar bias errors. The range 
bias error is the single most limiting error in the use of C-Band radar data for 
geodetic work. 

In the course of this study refinements made to the preprocessing 
procedures for conventional C-Band ranges significantly reduced their noise and 
systematic error. In particular, the AGC-based beacon delay variation is an 
important correction necessary in all future precision reduction of such data. 
Table 4.1 summarizes the calibration techniques developed or utilized in this 
study. 

The limiting factor in reliable elimination of systematic error in the 
integrated range data is the availability of very accurate estimates of the radar 
error parameters, especially the servo lag coefficient and ionospheric refraction 
effects. The elimination of all systematic error from a range data set with 
essentially negligible noise is a difficult task which will extend traditional radar 
calibration techniques far beyond the domain in which they are ordinarily 
applied. 
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TABLE 4.1. COHERENT C-BAND CALIBRATION TECHNIQUES 


NOISE REDUCTION 

• NON-RECURSIVE POLYNOMIAL FILTER TO REDUCE RADAR RECORDED DATA 
RATE (APPLIED TO RANGE AND RANGE RATE DATA). 

SYSTEMATIC ERROR REDUCTION 
RANGE DATA 

• PRE- AND POST-MISSION RANGING TO SURVEYED TARGET. 

• SIGNAL STRENGTH DEPENDANT BEACON DELAY VARIATION. 

• SERVO LAG CORRECTION. 

RATE DATA 

• TIMING ERROR (PRI/2). 

• ONE LINE UP/ONE LINE DOWN BIAS ESTIMATION TECHNIQUE. 

• INTEGRATION INTO PSEUDO-RANGE MEASUREMENTS. 
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APPENDIX A 


STEADY-STATE SERVO ERROR MODEL 


Radar tracking system output does not consist of instantaneous, 
uncorrelated measurements but rather the sampled output of a servo loop. 
Where a laser ranging device uses measured round trip time-of-f light as the basis 
for its output range measurement, the radar system uses the deviation of this 

quantity from its expected value (in the sense of range gate position) to update 

! ■ : - 

its servo. At any given instant in time the range servo is not likely to agree with 
the last instantaneous ranging "measurement," nor should it if its output is to 
have a reasonable noise level. 


All realizable (finite bandwidth) servomechanisms lag. The choice of 
bandwidth is usually a tradeoff involving noise reduction, dynamic response, and 
stability. Steady-state lag error in servomechanisms is generally modeled (I) as: 


6 : 6 \ 6 : 

«(t.) » T 1 + -! + ... 

K o K v K a 


where 0. and its derivatives represent the input signal or function to the 


servo. The constants K q , K v , and K q are servo design constants and figures of 


merit. Type 1 servos model only the input signal 9 , and lag to first order 
with 9. Type II servos, such as the DIRAM range machine and the Coherent 

i ' r—;-;- ' • • 

Signal Processor in the Wallops Island AN/FPQ-6, model both 9 and 9 , and !ag to 
• • 

first order with 9. Thus, the determining factor for first order lag in both 


systems is K q , as shown in Table A-l. (For both systems in all bandwidths K q is 
infinite and K y > I O'*.) (2) 








Figure A- 1 illustrates typical GEOS- 3 dynamics through a given horizon- 

to-horizon high elevation pass. The DIRAM range tracker lag error should go as 
! • • 1 
the R function* divided by K q , usually 200. In other words it will; lag 

approximately 25 cm at maximum acceleration. The Coherent Signal Processor 
operates on range rate and lags to first order as R (jerk) divided Jay its K , 
ordinarily 27. It would typically lag 1-2 cm/sec at peak R for a high elevation 
pass. The integrated CSP ranges, however, contain the integral of this In- 
dependent lag and hence, like the DIRAM ranges, will lag as R albeit with a 
considerably lower K Q of 27, corresponding to a 2 meter effect in a high 
elevation pass. 
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In conventional (O-C) residual space the pattern is reversed and the lag effect 
is "bucket shaped." 
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TABLE A— 1. 

FPQ-6 DYNAMIC LAG COEFFICIENTS 


DIRAM RANGE TRACKER 


BW SETTING 

1-5 6 

7 

8 

9 

10 

11 

12 

K a 

2.2 6.5 

20 

75 

200 

t 

1250 

1790 

3350 





1 

Normal 





Operation 


COHERENT SIGNAL PROCESSOR 

FINE LINE FILTER BW 5 Hz 15 Hz 40 Hz 160 Hz 

K_ 1 27 225 3000 

I 

Normal 

Operation 
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APPENDIX B 

LEAST SQUARES SOLUTION FOR ERROR MODEL TERMS 


The four parameter error model for integrated range residuals was fit to 
the residuals in the conventional least squares sense. The model 


4 R i = a o + a,R. + a 2 <, .-t 0 ) + a 3 R i 


where 


4R. = ith integrated range residual 


offset or bias; error in the constant of integration 


relative timing error between the CSP and DIRAM systems 


range rate measurement bias in the CSP system 


I 


R. 


1st order lag coefficient; ^ — ; see Appendix A 

a 

range rate at ith point 


R. 


range acceleration at ith point 


(tj-t Q ) = elapsed "integrator time" at ith point. 


leads to normal equations which in matrix form require: 


B- 


V* M . i, . i J I I I I I I • T 




The lower diagonal elements in the symetric matrix have been omitted 
for clarity. Computationally we find from the model 


34R, 


Sa 0 

34R. 


= I 


I 


= R. 


da 


34Rj 

3a 0 


(t. - v 


d4R. 


3 a, 


= R. 


The fitting process, implemented in an existing residual analysis program 
(GEORGE) of the GEODYN system, was thus largely a matter of summing the 
normal matrix entries and one 4x4 matrix inversion. The necessary range 
second derivative estimates (R) were produced by twice applying a fifth order 
central difference numerical interpolation scheme to the smooth, calculated 
range data available on the orbit determination program's output tape. 


